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Abstract 



"^J- | Starting from the first principles of nonrelativistic QED we have developed 

the quantum theory of the interaction of a two-component ultracold atomic 
ensemble with the electromagnetic field of vacuum and laser photons. The 
£> \ main attention has been paid to the consistent consideration of dynamical 

dipole-dipole interactions in the radiation field. Taking into account local- 
. field effects we have derived the system of Maxwell-Bloch equations. Optical 

t*"** ! properties of the two-component Bose gas are investigated. It is shown that 

the refractive index of the gas is given by the Maxwell-Garnett formula. All 
■ equations which are used up to now for the description of the behavior of 

an ultracold atomic ensemble in a radiation field can be obtained from our 
general system of equations in the low-density limit. Raman-Nath diffraction 
of the two-component atomic beam is investigated on the basis of our general 
^ . system of equations. 
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1 Introduction 



In recent years a great attention has been paid to the investigation of two-component 
Bose-Einstein condensates (BEC). A two- component BEC can consist of spatially 
separated identical atoms, or it can be a binary mixture of different alkali atoms, for 
instance, 87 Rb- 23 Na, or different isotops like 87 Rb- 85 Rb, or different hyperfine states 
of the same alkali atoms. A number of phenomena in two-component BECs, which 
are not possible in single-component BECs, has been theoretically predicted and 
some of them have been observed in experiments. It has been shown that the BEC in 
a double-well potential can oscillate between the wells by quantum coherent atomic 
tunneling [1], |2j. Oscillations of this kind can take place also in a two-component 
BEC, which consists of the same atoms in different internal states |§. Due to the 
nonlinearity arising from atom-atom interactions, the oscillations are expected to 
be supressed when the population difference of components exceeds a critical value 
in a process known as macroscopic quantum self-trapping (MQST) . However, in 
the process of collisions between the condensate and noncondensate atoms MQST 
decays away 0. The dynamics of spatial separation of two-component BEC has 
been studied in papers |5|, [(J . 

In the present paper we shall investigate optical properties of two-component 
BECs interacting with off-resonant laser radiation and develop mathematical for- 
malism for nonlinear atom optics with two-component condensates. Nonlinear atom 
optics with single-component condensates is a rather well studied subject. In pa- 
pers |7], H, |9|, [10, [ll[] different mathematical formalisms for the description of non- 
linear phenomena in atom optics of single-condensates has been proposed. Optical 
properties of the single-condensates subject to the influence of off-resonant laser ra- 
diation have been investigated in papers || [L_I], [TTj, |T2]. However, to our knowledge, 
nothing has been yet done in this direction for multicomponent condensates. Follow- 
ing the ideas, presented in our previous papers [[II], |ITJ , we shall derive the system of 
Maxwell-Bloch equations for nonlinear atom optics of two-component BECs. As an 
application of our general theory we shall consider a diffraction of two-component 
atomic beam from a standing light wave and discuss the specific features of this 
phenomenon, which does not take place in the analogous single-component process. 



2 The Hamiltonian for the two-component con- 
densate interacting with photons 

We consider a system of ultracold atoms which is a mixture of two species of two- 
level atoms with masses mi, m^, transition frequencies uui, u>2, and matrix elements 
of the transition dipoles moments d%, d 2 - We shall describe such a system in terms of 
matter field operators. Let \gj) and \ej), j = 1,2 are the vectors of the ground and 
excited states of the quantized atomic fields. Then the corresponding annihilation 
operators of the atoms in these internal states are ip g j and ip e j. Matter field operators 
are assumed to satisfy to the bosonic equal time commutation relations and the 
operators of different components are assumed to commute. 
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The Hamiltonian of the second quantized atomic field interacting with the pho- 
tons in the multipolar formulation of QED and in the electric dipole approximation 
can be written down in the following manner 

H = Ha + H F + Hai + Haf , (1) 
2 T r ( h 2 X7 2 



E J drflajfat) y-^^-J ^sj{r,t) + J dr$l j (r,t)fuj j 'ijj ej (T,t) 
H F = Y, huJ kci x (t)c kX (t) , Hm = - f drP(r,t)E in (r,t) , 

kA J 

H A f = - /drP(r,t)D mic (r,t) , 

where the operator of the microscopic displacement field is given by 



D mic (r, t) = 2^zd — r^e A c kA exp (tkr) + H.c. , (2) 
kA V v 

and the operator of the polarization field has the following form 

p = E Pi = E (p; + pj) = £ d i Wi4i + tiAi) ■ ( 3 ) 

i=i i=i 3=1 

Here we assume that the incident electric field E in (r, t) is produced by the laser, 
so it can be treated as a c-number function. In the Hamiltonian (|I|) we neglected 
all types of contact interaction. This approximation is valid when the saturation 
parameters of atomic transitions are small enough P, [TO]. We do not include into 



the Hamiltonian ([D a trapping potential, because our aim is to develop a theory of 
nonlinear atom optical processes of unconfined atomic beams. 

3 Heisenberg equations of motion for the atomic 
and photonic operators 

The Heisenberg equations of motion for the atomic and photonic operators are easily 
derived by from the Hamiltonian (JJ) and are given by: 

ift^M = _ d,E m (r,t)^,(r,t) (4) 

-ft E 9l\jcix(t) exp (-ikr) tp ej (r, t) - hip ej (r, t) E 9kXj exp (ikr) c kA (t) , 

kA kA 



dijj e (r,t) h 2 V 2 - 
ih — e ^ ' = _ ^~^i(r, t) + hujjip ej (r : t) - d i E in (r, t)i/) gj (r, t) (5) 



* ~t 

?kAj c k. 

kA kA 



~ h E ^kAjcLW ex P (-^kr) $ ffi (r, t) - fy£ OT -(r, t) E #kAj exp (ikr) c kX (t) 
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dc (t) 2 r 

ih m^ = huJkd ^ 9 *w J dre ^ r [^ (r ' t) ^ (r ' *) + ^ (r ' t) ^ (r > t} 

(6) 

where are the positive and negative frequency parts of the incident classical 
electric field. The operator products in Eqs.(f|),(|5|),([]) are taken in normally ordered 
form. 

The formal solution of (J|) for the photon operators is 

Ck\(t) = c kA (0) exp (-iu k t) + i V gl Xj [ dt' [ dr' exp [iu k (t' - t) - zkr'] (7) 

j=1 J° 

x [$.(r' ; t')^ gj (r', t') + ft gj (r', t')$ ej (r', t') , 

where the first term 6]<_\(0) refers to the free-space photon field and the second one 
goes back to the interaction with the atoms. 

To study the back reaction of the photons on matter we insert (|7J) in (^) and ([S]). 
By doing this procedure we eliminate photons in favor of atoms. In the rotating- 
wave approximation we obtain the following dynamical equations for the operators 
of the matter fields 

iti^^ = -^^ gj (v,t) - dB r Jv,t^ ej (r : t) : (8) 
d'i (r t) h 2 V 2 - 

in ito = ~l^T^ r ' ^ + h fa + ^' ~~ z ^' /2) ^' (r ' 

-^(r,t)dE+ c (r,t), (9) 

where 8j and jj are the Lamb shift and the spontaneous emission rate of a single 
atom in free space, respectively. We have introduced the operator of the local electric 
field 



EL( r > t) = E tn(*i t )+ i Y,\l ^y^ e AC k A(0) exp (ikr - iu k t) 



t /rt«v /^-"'» , (10) 

where Vx refers to the point r. The polarization operator P is given by eq.(§). Note 
that in Eq. ([370]) a small volume around the observation point r is excluded from the 
integration. 

Eq. (0) shows that E^ c (r, t) is a superposition of the incident field E in (r, t), 
vacuum fluctuations of the photon field, and the electric field radiated by all other 
atoms, which has exactly the same form as in classical optics. It is this local field 
which drives the inner atomic transition in Eqs.@, (§) which can be regarded as 
an atom-optical analogue of the optical Bloch equations [13, 14 1. They describe 
the dynamical evolution of second quantized matter in the field of electromagnetic 
radiation. 
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4 Lorentz-Lorenz relation and the system of Maxwell- 
Bloch equations in atom optics of two-component 
BEC 

4.1 Local-field correction 

The solution of Eqs. (|8]), (|9]) represents a rather complicated mathematical prob- 
lem because these equations contain explicitly dipole-dipole interactions. In many 
particular situations such a detailed microscopic description of matter is not nec- 
essary and it is more convenient to consider optical properties of the medium on a 
macroscopic level. This can be done by introducing the macroscopic field E mac (r, t), 
which satisfyes to the macroscopic Maxwell equations for a charge-free and current- 
free polarization medium, instead of the local field E; oc (r, t) in Eqs. (0). 
As in Ref. [|15[ we can introduce the macroscopic field by setting 



4-7T a 

E± (r,f) = E± ac (r,t) + y P ± (r,t) . (11) 

This equation is often called in the literature the Lorentz-Lorenz relation. It consti- 
tutes the basis of the local-field effects in classical |16| , quantum [|HJ and nonlinear 



optics (see [II, 18, [19] and references therein). In the case of a classical electro- 
magnetic field interacting with a macroscopic dielectric medium this relation can be 
derived from first principles under the assumption of homogeneity and isotropy of 
the dielectric medium. We take it here as the definition of E^ ac (r, t). It can then be 
shown with Eqs.(|3|) and (|T0| ) that this E^ ac (r,t) satisfies the macroscopic Maxwell 
equations, which can be written down in the form of the wave equation 

VxVxE^) = -^%^-f^. (12) 



c 



so it is justified to call it the quantum field operator of the macroscopic electric field. 
At the same time this definition allows us to interpret our results on ultracold atomic 
gases in analogy to the interaction of light with a macroscopic dielectric medium. 



4.2 Nonlinear matter equation 



We substitute ( |ITD in (||) and (||) and pass to a reference frame rotating with fre- 
quency ujl of the incident field, which is assumed to be monochromatic, to obtain 



ih 



dt 



2V72 



2mi 



h 



si 



- Y^10ll^ffl0el - ydida^L^/Wel , (13) 



ih- 



dt 



-h (Ai + 27l/2) (pel ~ — dld 2 '0 fl l'0g20e2 , 



(14) 
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in ~ftf = ~~2^^ 92 ~ ( r ^ e2 _ Y rf 20L^ 9 20e2 - ^-d 1 d 2 }li'}gl}e2 , (15) 

= - "2m2~^ e2 ~ 2 2 ^ - y«2^2^20e2 

47T - + - 

-H (A 2 + «72/2) 0e2 - —d 1 d 2 lpg2% 1 (f)el , (16) 

with the detunings Aj = u>l — ojj — 5j, j = 1,2. The position dependent Rabi 
frequencies ^f(r) = 2dj£^ ac (r) /% are related to the macroscopic electric field. 

Because we are mainly interested in atom optical problems and want to study 
the coherent evolution of the center-of-mass motion of the gas, we shall neglect 
spontaneous emission. This is valid for situations where the absolute values of the 
detunings are much bigger than the spontaneous emission rates and Rabi frequencies 
I A? I ^ 7? ' m or der to do this approximation self-consistently we drop in the 
following the vacuum fluctuations and the spontaneous emission rates jj from our 
equations. In addition we shall replace all the operators by macroscopic functions. 
We may, therefore, apply the adiabatic approximation ||, [13|, to fll4]), (|16D, which 
gives 

2A' oc (r,t) 



where the local detuning is given by 

47T 



Af (r,f) = Aj {l - ^ [ ttl |^i(r,t)| 2 + a 2 |^ 2 (r,t)| 2 ]} . (18) 

Here acj = —d 2 /hAj is the atomic polarizability for j-th component. 

Then substituting flT7D in (|T3|), (fH)|), which eliminates the excited states, we 
obtain as the result a system of nonlinear equations for the ground state matter 
fields ip g j(r, t) 



.^di> gj (r,t) h 2 V 2 hAj , 



A l ° c (r,t) 



Varying the parameters in Eq.([T^) we can change the nonlinear potential which 
is given by the second term on the r.h.s. of Eq . (|i~9|) . For instance, for increasing 
densities and positive detunings Aj the local detunings grow and, correspondingly, 
the nonlinear term in (|19|) representing the coupling to the macroscopic electric field 
becomes smaller. On the other hand, for negative detunings the absolute values of 
the local detunings decrease with increasing densities and the nonlinearity becomes 
greater. This behavior is exactly the same as we had in a one-component medium. 
In a two-component medium another regime is possible which can not be reached 
in a one-component medium: If the signs of the detunings are different, then one 
can increase the densities of the components in such a manner that the values of the 
local detunings, and therefore of the nonlinear potentials, will remain constant. 
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While Eq.([L9]) will allow us to derive an expression for the dielectric susceptibility 
of a Bose gas which closely resembles that of a classical gas we have to remark that 
it is only valid for low enough values of parameters Ej = atj \ip g j\ 2 - The reason is that 
the adiabatic approximation (|T7]) represents the first-order term in an expansion in 
1/Aj ||13|| . Therefore one also should expand Eq. ([19]) to first order in this parameter. 



This procedure leads to a pair of coupled Gross-Pitaevskii equations 



ihZ ryj_ = ,___ + 



dt { 2mj AAj 



~Y («1 \lpglf + a 2 \lpg2 2 



i/> gj . (20) 



Equations of this type have been used, for instance, in papers || f21|j . 

4.3 Optical properties of the two-component ultracold gas 

Making use of the adiabatic solutions ( |TTD we obtain the following expression for 
the medium polarization 

P+(r,t) = X (r,t)E+ ac (r,t), (21) 
where dielectric susceptibility is given by 

/ i \ _ J2 2 j=iaj\^gj(r,t)\ 2 

Dielectric susceptibility is a rather important parameter, because it describes the 
propagation of the laser radiation inside a medium. In most of the practical situa- 
tions the electromagnetic processes are much faster than the center-of-mass motion 
of the atoms. Therefore, \ can De considered as a time-independent quantity. Let 
us assume in addition that the spatial variations of the atomic density are not very 
large, such that Vx 0. Then divE^ ac « 0, and we have the following Helmholtz 
equation for the macroscopic electric field 

V 2 ^ ac + k 2 L n 2 £± ac = , (23) 

with the refractive index n given by the Maxwell-Garnett formula 

o i + teLi a, ivy 2 / x 

n 2 = l + An X = l^- 1 J \ W , (24) 

which is a two-component analogue of the Clausius-Mossotti formula. 

Eqs. dl9|) , (^), (p4]) can be considered as an atom optical analogue of the system 



of Maxwell-Bloch equations. In general they have to be solved in a self-consistent 
way and in usual situations solutions can be obtained only by doing numerical 
calculations. In the next section we shall consider one particular example of an 
analytical description of a nonlinear atom optical system. 
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5 Diffraction of a two-component ultracold atomic 
beam from a strong standing light wave 

We consider a typical scheme for the observation of diffraction in atom optics: An 
incident atomic beam moves in ^-direction, perpendicular to two laser waves counter 
propagating along the y-axis with wave vectors +ki and — k^, respectively, and 
with Gaussian envelope. From the uncertainty relation it follows that in order to 
get a distinct diffraction pattern, the width of the atomic wave packet w y should 
be sufficiently large compared to the wave lengh of the laser radiation in a medium. 
In this case the atoms can be described as a homogeneous medium with constant 
refractive index. If the spontaneous emission does not make any contribution, the 
effect of the atoms on the laser beam is purely dispersive and only the wavelength 
will be shifted. This means that in a medium we shall have a standing wave which is 
formed by counter propagating laser beams with the wave vectors +nk^ and — nk^, 
respectively. In this approximation the solution of (E3|) with fl2~4|) is given by 



|f2j| 2 exp z 2 /wf) cos 2 nk^y 



(25) 



We assume that the longitudinal kinetic energy of the atomic beam, associated 
with the center-of-mass motion in z direction, is large compared to the nonlinear 



potential in eq . fll9|) . Then the z-component of the atomic velocity will not change 
much and, therefore, the motion of atoms in z direction during the whole evolution 
can be treated classically. Only the motion in y direction should be treated quantum 
mechanically. In such a situation the coordinate z plays the role of time and we can 
change the variable t = z/v g j in ([TJ) with v g j being the group velocity of the j-th 
component. In addition we assume that we are in the Raman-Nath regime and 
we can neglect the transverse kinetic energy during the interaction of the atoms 
with the electromagnetic field. This approximation is valid for heavy atoms or if 
the interaction is so strong that atoms can take up momentum without changing 
considerably the velocity p2 |. In this case the density of atoms remains unaltered, 
but their phase changes. Making use of all these assumptions we can write down 
the solutions of eqs. (|T9|) for z 3> Wl (in the far zone) in the following form 



^gjiVi oo) = $gj(y, -oo) ex P 



V 



Vl(y, 



4A j v gj (l + V lPgl + V 2Pg2 y 



jdz 



(26) 



where 



Air 



-a. 



4tt og_ 
3h A~ 



and pgj = \ip g j\ is the density of atoms in the ground state. 
We represent p g j as Gaussian wave packets with width w y 

p gj = Pj exp (-y 2 /w 2 y ) . 



(27) 



(2f 
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Then we substitute Eqs.(^) and (|28|) into Eg. (p6|) and take into account that the 
width of the atomic wave packet must be much larger than the wavelenght of the 
laser radiation, i.e., w y 3> 27r/n/cz,. After integration we get the following result 

oo 

^ gj (y, oo) = ^ (y, -oo) £ e i2 ^(-i) %(r,) , (29) 

q= — oo 

which is represented here in the form of a Fourier series expansion. We use the 
notations: 

73 = 2 9j / (1 + Vfo + V 2 p 2 f , a = ^^V* • (30) 

J q is the g-th order Bessel function. 

From the solution ( p9|) it follows that the momentum transferred from the laser 
beam to the atomic beam is the same for both components and equals to 2qnki,y. 
It is determined by the wave number of the incident laser radiation k^ and the 
refractive index of the gas n. However, the probabilities to find the components of 
the beam in a momentum state 2qnkL are different for different components: 

P q] = J 2 q {r J ), g = 0,±l,±2,..., (31) 

with Poj being the probability to find the j-th component in the same momentum 
state as for the incident atomic beam. The angle of diffraction a q j for a particular 
momentum state q and for a particular component j is thereby given by 

2qnhki 

tana„j = . (32) 

rrijVgj 

Therefore the diffraction pattern, as it follows from Eqs.(p9|), (|31|), (|32"D, depends 
on the densities of the components. Depending on the values of mj and v g j the 
angle a.\ q can be either the same as a q2 or different. Only if m\V g \ = m 2 v g 2, i.e., 
when the momenta of different components associated with the group velocities are 
the same, a q i = a q2 . In all other situations, for instance, if the group velocities 
of the components are equal to each other or if we have a monoenergetic atomic 
beam, a q \ 7^ a q2 , and in the diffraction pattern one can observe spatially separated 
components. 



6 Conclusion 

In the present paper we have investigated the process of the interaction of a two- 
component BEC with the field of vacuum and laser photons. The two-component 
BEC is treated as a binary mixture of two-level atoms with different masses, tran- 
sition frequencies and transition dipole moments. Starting from the microscopic 
model and making use of the multipolar formulation of QED, a general system of 
Maxwell-Bloch equations is derived which can be used for the description of nonlin- 
ear phenomena in atom optics. Optical properties of the two-component BEC are 
investigated. The refractive index is shown to satisfy the Maxwell-Garnett formula. 



9 



As a typical atom optical application, we have considered the diffraction of two- 
component atomic beam from a strong standing laser wave in the Raman-Nath 
approximation, which allows to obtain simple analytical solutions. It is shown that 
in most of the situations one can observe splitted components of the beam in the 
diffraction pattern. 

The limits of validity of the results, obtained in the present paper, are essentially 
restricted by the adiabatic approximation, which is correct up to the first order with 
respect to the small parameters 1/Aj. Therefore, our results are valid for small 
enough Ej = a>j \ip g j\ 2 ■ They generalize our previous results JTD|, |TT |. 



Although, we have considered here explicitly only a two-component BEC, the 
generalization to an arbitrary number of different atomic species is straightforward 
and can be done very easily. 
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